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The hot nuclear matter created at the Relativistic Heavy Ion Colhder (RHIC) has been char- 
acterized by near-perfect fluid behavior. We demonstrate that this stands in contradiction to the 
identification of QCD quasi-particles with the thermodynamic degrees of freedom in the early (fluid) 



stage of heavy ion collisions. The empirical observation of constituent quark 



scaling of elliptic 



flow [l[ is juxtaposed with the lack of such scaling behavior in hydrodynamic fluid calculations fol- 
lowed by Cooper-Frye freeze-out to hadrons. A "quasi-particle transport" time stage after viscous 
effects break down the hydrodynamic fluid stage, but prior to hadronization, is proposed to rec- 
oncile these apparent contradictions. However, without a detailed understanding of the transitions 
between these stages, the "n," scaling is not a necessary consequence of this prescription. Also, if 
the duration of this stage is too short, it may not support well defined quasi-particles. By comparing 
and contrasting the coalescence of quarks into hadrons with the similar process of producing light 
nuclei from nucleons, it is shown that the observation of "n," scaling in the final state does not 
necessarily imply that the constituent degrees of freedom were the relevant ones in the initial state. 

PACS numbers: 25.75.Dw 



INTRODUCTION 

The understanding that hadrons are not fundamen- 
tal particles, but rather composite systems comprised of 
asymptotically- free quarks and gluons, led directly to the 
idea that hadronic matter under extreme conditions of 
temperature and pressure would transform into a "quark- 
gluon plasma" In this state, quarks and gluons could 
be thought of as quasi-free particles moving along trajec- 
tories, such that the thermodynamic properties could be 
approximated as that of an ideal gas with non-interacting 
degrees of freedom. More specifically, quarks and gluons 
would no longer be confined in their parent hadron, but 
would be free to roam about larger volumes. This sce- 
nario gives 40 degrees of freedom for the quark and gluon 
states (including two light quarks [HI]) as opposed to 3 
pion degrees of freedom in a hadronic system at a lower 
temperature. This enumeration is often associated with 
the order of magnitude change in the entropy s/T^ and 
energy density e/T^ as observed near the transition tem- 
perature in lattice QCD calculations |3i]. 

This physical picture has generally been thought to be 
naive, since the high densities and moderate momentum 
transfers should lead to short mean free paths. However, 
it has a certain persistence in the heavy ion literature 
(e.g. the CERN press release [1] or the discussion of such 
in the review as it presents an easily-described phys- 
ical scenario with nominally well-defined consequences. 
Over the last few years, experimental data has suggested 
a paradigm shift where the quark-gluon plasma is de- 
scribed as a nearly inviscid fluid, and yet the picture 
of quasi-particles carrying the thermodynamic degrees 



of freedom often persists. Notably, recent experimental 
papers describing the scaling properties of elliptic flow 
measurements state that "scaling suggests that partonic 
collectivity dominates the transverse expansion dynam- 
ics 6]" and "[scaling indicates] a pre-hadronization state 
in which the flowing medium reflects quark degrees of 
freedom Q." The analysis presented in this paper sug- 
gest that these conclusions are not well supported either 
by experimental data or theoretical ideas. 



QUASI-PARTICLES IN THE QUARK-GLUON 
PLASMA 



Interactions between propagating quarks and gluons 
in the quark-gluon plasma generates thermal masses, 
and thus the quarks and gluons can be called "quasi- 
particles." A quasi-particle's mass is determined by 
the local properties of the medium, which thus distin- 
guishes between Lorentz frames. By contrast a "real" 
particle has a Lorentz-invariant rest mass. Such quasi- 
particle descriptions allow even strongly coupled systems 
to show properties reminiscent of free non-interacting 
systems. The characteristic feature which indicates so- 
called "good" quasi-particle states is a decay width 
(which characterizes the coupling strength and number 
of interaction channels) that is smaller than its mass. 

The dynamical role of quasi-particles in the hot phase 
of QCD is under intense debate. In [sj an accurate match 
to lattice QCD calculations of the entropy for a purely 
gluonic medium is obtained in a perturbative calculation 
involving quasi-particles for temperatures down to 2Tc. 
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They note that "ahhough the quasiparticle picture sug- 
gested by such fits is a rather crude representation of the 
actual physics of non-abehan gauge theories, it supports 
the idea that one should be able to give accurate descrip- 
tion of thermodynamics of the QCD plasma in terms of 
its elementary excitations." Such quasi-particle models 
have been extended and applied down to l.OSTc where 
they are able to describe deviations from the free quark 
and gluon picture d, H^. However, they find that these 
effective degrees of freedom are rather heavy and have a 
sizable width making it difficult to associate them with 
"strict" quasi-particles. Thus, it is not obvious that the 
identification of the hot dense system with a quark-gluon 
plasma necessarily implies a quasi-particle picture, even 
where perturbative calculations are possible. 

Additional information on quasi-particles near the 
transition temperature are available from lattice QCD. 
The suggestion that the bulk thermodynamics of the 
quark-gluon plasma are governed by a large set of colored 
bound states of light quarks flll| has largely been ruled 
out by examining fluctuations of quark number and elec- 
tric charge [l2|. Additionally, baryon-strangeness corre- 
lations place severe limits on qq bound states [i3,[i3. It 
is notable that the consistency of Xud with a weakly in- 
teracting plasma of quarks and anti-quarks above Tc 14 1 
does not provide compelling evidence that these are well 
defined quasi-particles dominating the bulk thermody- 
namics. 

In many systems the thermodynamic degrees of free- 
dom can be counted by associating them with particles or 
quasi-particles and accounting for their thermal motion 
and quantum numbers (for example \kT and |A:T for 
non relativistic monatomic gases and diatomic gases re- 
spectively). A key question is whether or not this is the 
case for the quark-gluon plasma. The thermodynamic 
properties of some systems, such as Landau's Fermi liq- 
uid [3], are dominated by the quasi-particle degrees of 
freedom, while in other systems quasi-particles are a sin- 
gle medium excitation (for example the fractional quan- 
tum hall effect [Hj). In the case of hot QCD matter 
both may be interesting, but the former is the more rele- 
vant when discussing the quark-hadron phase transition 
observed on the lattice. 



FREEZE-OUT FROM HYDRODYNAMICS 

One of the great discoveries of the RHIC program has 
been that simple hydrodynamic models provide a good 
description of RHIC data, es pec ially the development of 
radial and elliptic flow, see l3| and references therein. 
The hydrodynamic equations are based on the energy- 
momentum tensor, and viscosity is included as an ex- 
pansion parameter in terms of velocity gradients. The 
viscosity term is set to zero for a mathematically per- 
fect fluid. When invoking a picture of quasi-particles, a 



connection can be made between the mean free path of 
these objects via kinetic theory and the viscosity term in 
the hydrodynamic equations. The hydrodynamic limit 
of the quasi-particle picture is reached when the typical 
gradients in the flow involve length scales that are large 
compared to the typical mean free paths. 

First-order estimates suggest that experimental data is 
consistent with essentially zero viscosity fis'] , but cannot 
rule out a viscosity to entropy density ratio consistent 
with the lower bound derived from the AdS/CFT dual- 
ity {rj/s > I/Att) [l^. Furthermore, estimates of collec- 
tive behavior based on classical transport (Boltzmann) 
approaches are generally unable to describe RHIC data, 
whether due to the small cross sections implied by a per- 
turbative transport picture [20], or the large formation 
times required by hadronic transport .21|. 

The hydrodynamic behavior of the medium must even- 
tually break down as the density of the medium drops 
below a threshold value, and particles "freeze out" into 
vacuum. In addition, there may be particles (for exam- 
ple at high transverse momentum) that are not equili- 
brated into the bulk medium during the limited evolu- 
tion time (i < 10 — 15 fm/c) ^2^. Thus, in order to 
compare hydrodynamic calculations with experimental 
data (i.e. measured hadron distributions after freeze- 
out), hydrodynamic calculations are typically terminated 
at a chosen temperature T . This is usually chosen to 
be the Hagcdorn temperature Tch ^ 170 MeV in sin- 
gle freeze-out models, and Tth ^ 100 MeV when a ther- 
mal freeze-out temperature is postulated due to hadronic 
re-scattering. In either case, this prescription deflnes a 
"freeze-out hyper-surface" in space-time, also referred to 
as a "surface of last scattering." The Cooper- Frye formal- 
ism [2^ is then used to decay fluid elements to hadrons 
by assuming local statistical hadronization. Many cal- 
culations also include the important non-equilibrium re- 
scattering between hadrons after hadronization up to the 
point of freeze-out [3, ll^ . 

Shown in Figure [1] are the published experimental 
data P, 0, 27 1 on elliptic flow V2 as a function of trans- 
verse momentum for various hadron species in minimum 
bias Au -f Au collisions at ^/sWn — 200 GeV. Also shown 
are the results of a particular hydrodynamic calculation 
terminated with Cooper- Frye freeze-out [28|. This cal- 
culation begins with an equation of state that assumes 
chemical equilibrium until kinetic freeze-out (EoS Q), 
and a first order phase transition at 165 MeV. The 
initial time is taken as t—0.6 fm/c, and freeze-out oc- 
curs when the temperature drops to r/o=130 MeV. The 
particle spectra and elliptic flow (U2) qualitatively agree 
with the experimental data for values oipx < 1.5 GeV/c. 

The reasonable agreement of these calculations with 
the data, even without an extensive parameter sensitiv- 
ity or error analysis, has generally been taken as evidence 
that collisions of nuclei at RHIC form a perfect invis- 
cid fluid (i.e. small rj/s) with short thermalization times 
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FIG. 1: (color online) Elliptic flow «2 for identified hadrons 
as a function of pr in minimum bias Au + Au collisions at 
y/SNN = 200 GeV, compared with hydrodynamic calculations 
assuming Cooper-Frye freeze-out. 



(tq < 0.6 fm/c). Notably, different calculations with dif- 
ferent equations of state or initial conditions do not all 
agree with the data at the same level Similarly 
small values of rj/ s have also been derived via measure- 
ments of charm suppression and flow fs^. Of course, the 
quantitative importance of viscosity has only started to 
get serious treatment, both experimentally and theoreti- 
cally [Ulliilll. 

It is common in the field for the full range of trans- 
verse momenta to be sub-divided into three categories: 
The "low pt" region [px < 1.5 GeV) where hydrody- 
namics agrees with experimental data; the "intermedi- 
ate pt" region {px ~ 1.5 — 4.5 GeV/c) where hydro- 
dynamics over-predicts elliptic flow V2', and the "high 
Pt" region (px > 4.5 GeV) where hadronization from 
jet fragmentation is thought to be dominant. In the "in- 
termediate pt" region, the violation of "hydrodynamic 
scaling" and the observation of enhanced (anti) baryon 
relative to meson yields (often referred to as the "baryon 
anomaly" [iiilil) led some to postulate a different 
hadronization mechanism. 



RECOMBINATION AND "n," SCALING 

It has recently been speculated that in RHIC colli- 
sions [sHI a thermal distribution of quarks might coalesce 
(a process sometimes referred to as "recombination") into 
hadrons when they satisfy the minimum quantum num- 
ber requirements. If these requirements are met in a 
local region of phase space, then only the valence quarks 
of the correct flavors (and no additional constraints) are 
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FIG. 2: (color online) The same data as in Figure 1, plotted 
as V2 vs. rriT — m, the hadron transverse kinetic energy, both 
scaled by the number of constituent quarks per hadron. The 
reference line is a fit to the kaon data. 



needed to form a hadron. In [33|, they argue that in 
the light-cone gauge the hadron wave-function can be ex- 
panded in a Fock space of quarks and anti-quarks where 
the first order term is simply the valence quark content 
(e.g. \p) = \uud) + \uudd(I) -I- . . .), and that this is appli- 
cable when pt ^ m for the hadron. This simple picture 
has been successful in qualitatively understanding the 
"baryon anomaly" and flow patterns. Interestingly, this 
"recombination" picture appears to involve a decrease 
in entropy; although, for "intermediate p^" hadrons one 
can always posit that compensating additional entropy 
is generated in the bulk medium at lower px- In fact, it 
has been suggested in [s^ that the lattice QCD equation 
of state permits isentropic hadronization. 

It has been empirically observed [ij that all hadrons 
follow a universal trend when V2 versus px is repre- 
sented instead as V2/nq versus {mx — m)/nq, where nq is 
the number of valence quarks in the hadron formed and 
mx = \/pfT 



m? is the transverse mass. This was shown 
with STAR and PHENIX data [l| (as reproduced here) 
in Figure [21 and expressed as a function of the transverse 
kinetic energy KEx — tot — "m. This scaling is often 
termed "constituent quark scaling." 

If there are quasi-particles with the quantum numbers 
of quarks that are coalescing into hadrons, then an im- 
mediate question is what are the values for the thermal 
masses and widths. The "constituent quark" nomen- 
clature implies an effective mass of mpf/i 300 MeV. 
However, the 300 MeV constituent quark mass in a pro- 
ton is generated by the spontaneous breaking of chiral 
symmetry {{qq) ^ 0). Lattice QCD results indicate an 
approximate restoration of chiral symmetry in the quark- 
gluon plasma region above the critical temperature [3^ . 
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FIG. 3: (color online) The same data as in Figure 2, divided 
by a reference function fit as described in the text. 



FIG. 4: (color online) The same hydrodynamic calculations 
shown in Fig. 1, shown as V2/nq vs. {rriT 
absence of the scaling seen in the data. 



m)/nq. Note the 



Any such quasi-particles would be expected to have an 
effective mass, but not from chiral symmetry breaking. 
Therefore the hght quark masses assigned in [38| seem a 
remarkable coincidence since the dynamical mass of order 
gT that might be generated in the thermal medium does 
not violate chiral symmetry and is quite different from 
a constituent quark mass. In addition, since the scaling 
appears to work for mesons (including the pions, which 
are anomalously light Goldstone bosons), a simple coa- 
lescence of a 300 MeV quark and a 300 MeV anti-quark 
easily violates local energy-momentum conservation. 

If the quasi-particles containing different flavor quan- 
tum numbers (e.g. light quark flavor versus strange) 
have different masses, some scaling violations are ex- 
pected [3^. However, in the limit that the mass of such 
quasi-particles is much less than their momentum, the 
scaling with Uq gives no information about the mass dif- 
ferences. In the case of charm flavor quasi-particles, the 
current quark mass difference (~ 1.4 GeV) is always large 
compared with gT and the momentum up to intermedi- 
ate pt, so future measurements of scaling agreement or 
violation are particularly interesting 3^, 40| . 

Despite the categorization of the dominant physics in 
different pT ranges, the scaling does seem to work well 
at all PT- To examine this quantitatively, the data is fit 
by a reference curve of the form vijnq = atanh(/3(TOT — 
m)/nq) (as shown in Figure [2]), where the parameters a 
and [3 are determined from the charged kaon data alone. 
Fitting the kaon data is arbitrary and thus when dividing 
by the reference line, one only obtains a measure of the 
magnitude of deviation of the different hadrons from each 
other as shown in Figure O 

The scaling behavior for (tot — to) /nq appears to hold 
at the ±5% level for all hadron species measured for 
(tot - rn)/nq > 0.2 GeV. Note that (tot - ■m)/nq = 0.2 
GeV corresponds to pT — 0.52, 0.75, and 1.21 GeV/c for 
pions, kaons and protons respectively. Below this value. 



the deviations are larger (e.g. ± %15"25%); It is notable 
that the charged and neutral kaons have the largest devi- 
ation from each other, which may indicate some experi- 
mental issues. We note that there are first measurements 
of other hadrons (for example 0, 0| and and light 
nuclei d and '^He [g, |4l[), but the current errors only al- 
low one to confirm general agreement with the scaling 
behavior. 

The "rig" scaling at low px is not a natural prediction 
of hydrodynamics followed by Cooper- Frye freeze-out. 
These calculations depend only on the particle masses, 
and not on the quark content. However, it has been noted 
that if there is a regime where V2 is linear as a function of 
Tot — TO, then any re-scaling V2/a and (tot — m)/a will 
work for any constant a [42]. In order to test this, in Fig- 
ureUl a hydrodynamic calculation is plotted as a function 
of the scaling variables vijuq and (tot — m)/nq. Again, 
the different hadron species have been compared to the 
kaons, representing an arbitrary choice of reference. It 
is notable that the hydrodynamic calculation results are 
sensitive to the choice of initial conditions and equation 
of state. However, we do not expect these variations or 
the addition of hadronic re-scattering to modify the scal- 
ing results. The results reveal ±15-25% level deviations 
of the various hadrons from the reference line; whereas, 
the experimental data shown on the identical scale in 
Figure [3] do not. The hydrodynamic calculations do not 
have a linear relation between these two variables even 
at low pt, as a function of either tot — to or pt, at a 
level that precise experimental data has tested. 

As seen in the Figure [31 it appears that the data follow 
the Uq scaling hypothesis within experimental uncertain- 
ties above (tot — m)/nq « 0.2 GeV. The hydrodynamic 
calculation shown in Figure[4]does not, particularly in the 
so-called "hydrodynamic regime." This observation does 
not necessarily imply that a full recombination calcula- 



tion would follow the "riq" scaling as well as the data, 
but it is clear that hydrodynamics fails to predict this 
emperical observation. 

RESONANCE CONTRIBUTIONS 

The level of agreement of the elliptic flow data with 
"n^" scaling is quite surprising, considering all of the rea- 
sons that exist for it not to work (for example see 
In fact, if all hadrons are formed via recombination of 
constituent quarks, then one might expect that p mesons 
and A resonances would follow the scaling, but their de- 
cay daughter products might not S| . The decay of these 
resonances results in pions with a distorted V2 scaling, 
since the decay blurs the emission angle and shifts the 
Pt of the pion relative to that of the parent resonance. 
We have quantified this effect using a Monte Carlo sim- 
ulation. We show in Figure [5] the vijuq for tt from p 
decay, tt from A decay, and pions from all sources com- 
bined. In the simulation. A, p and tt particles are given 
a -pT according to the spectra calculated in [i^ with hy- 
drodynamics. In our simulation, the p and A resonances, 
and the directly hadronizcd non-decay pions are assumed 
to follow the "universal" scaling curve (reference line) 
drawn in Figure [2l The p and A resonances are sub- 
sequently allowed to decay, and the flow of the pions is 
computed. Shown in Figure [6] are the simulation results 
divided by the original reference line. Interestingly, the tt 
from A decay has a significantly larger V2 than the scal- 
ing, as previously noted in 44]. In fact, by construction, 
the elliptic fiow for tt from A decay must asymptotically 
approach 1.5 times the thermal pion value at large pT 
where the decay products are essentially co-linear with 
the parent. This large deviation in the "rig" scaling is de- 
emphasized in the combined result (iXaii ) due to the small 
fraction of A decay pions in the sample. In contrast, we 
find that tt from p decay have a slightly reduced be- 
low (m,T — m)/nq < 0.7 GeV. In fact, when the different 
contributions are combined, the scaling law appears "ac- 
cidentally" obeyed at the 10% level with a deviation in 
a region of (m^ — m) /uq similar to that observed in the 
experimental data. 

QUASI-PARTICLE IMPLICATIONS 

Assuming that quasi-particles with the quantum num- 
bers of quarks exist in the bulk presents some theoret- 
ical interpretation challenges since the light-cone gauge 
calculation assumes hadrons with px ^ m and moves 
the entropy issue to lower px- In addition, it has been 
noted by many authors that the scaling is not an obvious 
consequence of hadronization via quark coalescence, but 
instead a consequence of the phase space distributions of 
the recombining quarks, e.g. as shown in [43]. 
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FIG. 5: (color online) Simulation results for pion vz from 
different sources are shown. 
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FIG. 6: (color online) The results from Figure 5 are divided 
by the reference line. 

One might naively conclude that the medium is a per- 
fect fluid and is composed of flowing quasi-particles with 
the quantum numbers of quarks (though yet unidentified 
mass and width characteristics). Perhaps these quasi- 
particles even carry the thermodynamic degrees of free- 
dom in their particle form. However, this picture does 
not sit well with the hydrodynamic paradigm of early 
time dynamics at RHIC. In the case of a perfect fluid, 
there can be no well-defined quasi-particles, as their fi- 
nite mean free path would induce dissipation and viscous 
effects. 

Alternatively, if the medium is not described by a 
mathematically perfect fiuid (i.e. rj/s — 0), but instead 
one with a viscosity to entropy density ratio very near 
the AdS/CFT bound, -q/s l/Air, one might wonder 
whether or not it is then possible for the medium to con- 
tain well defined quasi-particles. Although a rigorous 
answer to this question is in principle difficult to find, 
it is possible to derive order of magnitude estimates for 
the ratio of a quasi-particle's width F, to it's mass, m 
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from basic relations. As mentioned previously, a "well- 
defined" quasi-particle is a mode that can be character- 
ized by a width much smaller than its mass, F/m ^ f. 
In the case of a non-relativistic dilute gas [45|, one finds 
that the viscosity 77 is proportional to the product of the 
average particle momentum, (p) = mv, the number den- 
sity, n, and the particle's mean free path. A: 



1 



V 



(1) 



In analogy with the thermodynamics of a non-interacting 
gas of massless bosons, the entropy density is assumed to 
be proportional to the number density, as s ^ An. 

To introduce the width of the particle, we identify the 
mean free path with the particle's lifetime, r, so that 
A ~ VT, and note that F and r are related by F = 2/r. 
Taking the ratio of F to m with these relationships in 
place suggests that the quasi particle width is propor- 
tional to the temperature, and similar to its mass: 



F 

m 



16T 
3m 



(2) 



where we have used the Maxwell velocity v = \/%T jixm 
to characterize the average quasi-particle speed. 

An alternative estimate begins from the expression for 
the viscosity to entropy density ratio for a weakly coupled 
quark-gluon plasma derived from kinetic theory in [46|, 



(!) 



SJ wQGP 



XT 



(3) 



Using again the relationships between A, r, and F, it 
happens that the width to mass ratio is once more on 
the order of the temperature, this time 



F 

m 



SttT 
15m 



(4) 



where we have assumed in this case that v ^ c/3 — 
In both of these estimates, if we assume that the thermal 
mass of quasi-particles is 1 — 3 • T, then the ratio F/m 
is of order one. Thus, if the produced medium has a 
viscosity near the conjectured bound, it is unlikely that 
quasi-particle modes would be well defined. 

Recently, it has been speculated that the quark-gluon 
plasma might have a small viscosity and maintain well de- 
fined quasi-particles, if there is a dominant contribution 
from turbulent color fields to the transport coefficients, 
referred to as an "anomalous viscosity" [47, 48']. This 
scenario allows for weakly coupled quasi-particles (and 
thus well-defined) and a large coUisional viscosity Typ, but 
where the total viscosity is determined by the anomalous 
term, i.e. 77^^ = ry^^ + 77^^. We note however that the 
QED plasma analog is typically non-relativistic, and in 
the QCD case of the quark-gluon plasma the system is 
relativistic. One could imagine that the heavy flavor el- 
ements of such a QGP plasma would be non-relativistic 



and thus the compatibility with the QED plasma calcu- 
lations would be more obvious. 

In a relativistic quantum field theory, the sharp dis- 
tinction between fields and particles does not exist and 
the factorization of viscosity terms above is non-trivial 
to mathematically define. Thus, while it is very interest- 
ing to experimentally determine the possible size of this 
anomalous viscosity term, it is unclear if it reconciles 
the perfect fluid and quasi-particle pictures. Applying 
non-relativistic calculational insights to the relativistic 
quark-gluon plasma, including the previous derivations 
for F/m, is common. However, it is often unclear how to 
determine which physics effects are actually relevant for 
describing the observed data. 

While quasi-particle widths are large (e.g. while T > 
l.OSTc, which happens to be close to the Tch assumed 
in most calculations), then hydrodynamics should be the 
appropriate description of the collision evolution, as is 
already assumed for RHIC collisions. However, as the 
temperature approaches Tc from above, the widths may 
become small enough at a temperature T^p, such that 
below this, a quasi-particle transport (QPT) approach 
would be applicable, with quasi-particles following clas- 
sical paths. The transition from fluid to quasi-particles 
would occur at the T — Tqp hyper-surface with the fluid 
V2 at that moment. At the end of this this QPT stage, the 
quasi-particles could form hadrons via recombination. 

This quasi-particle transport scenario preceding 
hadronic transport and then freeze-out allows the simul- 
taneous use of hydrodynamics (at T > Tqp), thermal 
hadron formation (for Tch < T < Tqp) and even hadronic 
re-scattering (for Tth < T < Tch, e.g. as implemented 
in Ref.[24]) while allowing in principle for "rig" scaling. 
With better knowledge of the quasi-particle masses and 
quantum numbers it could be possible to implement a 
multi-stage dynamical model incorporating all of these 
with realistic transport. Of course, the scaling is not 
inevitable, but depends quite sensitively on the details 
of the dynamical evolution. For example, if QPT takes 
place too early, and freeze-out to hadrons occurs rapidly, 
then hadronic re-scattering may well destroy the "n^" 
scaling that could have been present just after hadron 
formation. In addition, if the QPT stage is too short, 
the uncertainty principle would in fact preclude narrow 
width quasi-particles. 



HADRONIC TRANSPORT AND 



SCALING 



In '49'!, the string and hadron cascade models RQMD 
and UrQMD provide another means to understand "riq" 
scaling without hydrodynamics or degrees of freedom car- 
rying valence quark numbers. They consider a picture 
far from the inviscid fluid limit, comprised of hadrons 
with long formation times (r w 1 fm/c) and transport 
dynamics driven by various in vacuo interaction cross 
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sections. The supposition is that hadron re-scattering 
encodes the scahng with valence quark number through 
different Additive Quark Model (AQM) hadron-hadron 
cross sections. The use of AQM for hadronic cross sec- 
tions in these transport models is directly motivated by 
the concept that each constituent quark carries a fraction 
of the total hadronic cross section. Thus, the constituent 
quark number scaling it built into these cross sections. 
It is then immediately seen that "n^" scaling does not 
require the relevant degrees of freedom to be carried by 
separate quark-like objects. 

Two significant issues pertain to this picture. First, 
the overall magnitude of the V2 is substantially under- 
predicted by these calculations. Of course, if the forma- 
tion time is reduced, the final V2 value can be made to ap- 
proach experimental data. As one approaches the hydro- 
dynamic limit where the mean free path is much shorter 
than the system size, the relative size of the cross sec- 
tions no longer contribute significant differences in flow. 
In the hydrodynamic limit, all cross sections need not be 
the same, but rather, the corresponding mean free paths 
must all be less than some limiting value. Additionally, 
in [4^ they excluded the (j) and fl particles, which have 
OZI suppressed hadronic cross sections. In this case, the 
expected cross section is quite different in the scenario 
of separate quark-like objects from that where they are 
bound together. Thus, these particles actually represent 
an excellent test of their hypothesis. Another additional 
test would be to show the flow for heavy flavor particles 
(such as D mesons) jsof . 

We noted earlier that if one shortens the hadron for- 
mation time, one can reproduce the overall magnitude 
of elliptic flow. If such a calculation were to reproduce 
the experimental data for all hadrons, it might be tempt- 
ing to conclude the medium dynamics are dominated by 
hadronic transport, or at least that an ambiguity ex- 
ists that precludes ruling out this scenario. However, 
in- medium broadening of resonances is likely to occur, 
which is due to short time scale interactions (so called 
"collision broadening"), see for example [lOl and refer- 
ences therein. In the case of multiple interactions of res- 
onances on the time scale of 0.1 fm/c, the widths would 
be comparable to the masses and thus the calculation 
would not be self-consistent in it's treatment of hadrons 
in a cascade approach. 



LIGHT NUCLEI FLOW 

We have also performed a Monte Carlo coalescence cal- 
culation for light nuclei and their resulting elliptic flow. 
In contrast to quark recombination, the mechanism of 
nucleosynthesis in heavy ion reactions is well understood 
and the light nuclei wave- functions are known. We pa- 
rameterize an exploding fireball freeze-out by a Lorentz 
boosted distribution of nucleons with Boltzmann momen- 
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FIG. 7: (color online) Coalescence calculation comparing 
V2/nA for nucleons (N) and deuterons (d) and '^He for mini- 
mum bias Au + Au collisions at ^snn = 200 GeV. Prelimi- 
nary STAR data is shown for d (squares) and ^He (triangles) 
and PHENIX data is shown for (anti) protons (circles). ,41] . 



tum and an induced anisotropic flow 
these nucleons as in 



We then coalesce 



5l[ , and compute the flow of light 
nuclei shown in Figure [71 The experimental data for 
(anti) protons do not follow the reference curve perfectly. 
This is because the reference line is from Figure 2 and is 
fit to the charged kaons. 

The calculation results for deuterons ('^He) differ from 
the simple scaling of V2/A versus (m^ — m)/A of order 
10% (20%) at high uit — rn. In fact, this is expected as 
can be analytically shown in the case of a delta-function 
in momentum space as the coalescence condition. 



2< 



(pt/2) 



2v^{pt/2) (5) 



1 



2< 



{pt/2Y 



This is the same as the result from analytic quark coales- 
cence calculations [s^ |4o!] . Given the large V2 for protons 
at large tot — to, the correction term in the denominator 
exactly accounts for the deviation between the calculated 
curves shown in Figure [Tj 

In the case of light nuclei coalescence where the wave- 
function and mechanism are relatively well understood, A 
(or equivalently riq) scaling appears to be reasonably well 
obeyed simply by recombining nucleons close in space- 
time and momentum. In this example it would be incor- 
rect to say that the observed scaling (for example of light 
nuclei V2/nq with (mT—m)/nq or V2/A with (mT—m) /A) 
stems from baryons being the objects flowing with the 
medium when the flow is built up. This is not a logical 
conclusion from the scaling. Rather, the baryons happen 
to be flowing with the medium at a particular time, and 
the coalescence mechanism imprints that pattern on the 
light nuclei. 

The situation is quite parallel to quark coalescence in 
that one cannot conclude that constituent quark degrees 
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of freedom were relevant as the flow developed. One only 
needs these quasi-particles to have the relevant valence 
quark degrees of freedom and to have the right azimuthal 
anisotropy as the fluid breaks up. An interesting ques- 
tion is whether one can at minimum conclude that light 
nuclei (in the light nuclei coalescence case) were not the 
degrees of freedom at earlier stages. This is not obvi- 
ous. Just because the coalescence mechanism is consis- 
tent with the observed phenomena does not prove that 
coalescence is the source of light nuclei. For this con- 
clusion to be valid one would have to show that all the 
features present in light nuclei production are consistent 
with the coalescence mechanism and that other models 
for nucleosynthesis fail to describe the empirical results. 
Similarly to the quark coalescence case, can we exclude 
bound states (whether hadrons or something more ex- 
otic) as the relevant degrees of freedom at earlier stages? 
The previously discussed UrQMD calculations and un- 
certainties in the hadronization mechanism leave this as 
an open question. 

Light nuclei were formed in the early universe by pro- 
cesses such as p + 71 — > c? + 7, where the 7 is required to 
conserve energy and momentum. In heavy ion reactions 
these processes have cross sections too small to account 
for the yield of light nuclei [51]. In this case, the coales- 
cence prescription assumes only that the p and n must 
overlap with the deuteron wavefunction, and the remain- 
ing exchange of energy and/or momentum with the sur- 
rounding medium has little impact. In fact, deviations 
from a common penalty factor for each additional coa- 
lesced nucleon has been found to scale with the binding 
energy per nucleon in the nucleus [H^ l- This may simply 
reflect the larger capture cross section for more tightly 
bound light nuclei, or may be indicative of the required 
energy and/or momentum exchange with the medium. 
In the case of quark coalescence, the possible required 
exchange of energy and/or momentum with the medium 
is much larger. If the quark-like quasi-particles q have for 
example mass 1.5T w 255 MeV, then the formation pro- 
cess q + q — > TT requires a large transfer of energy and/or 
momentum to the medium. In fact, for all hadrons the 
exchange is much larger than in the case of light nu- 
clei. At the transition temperature, the mechanism of 
this exchange as the chiral condensate is re-appearing is 
not known. As pointed out in 42], the scaling of V2/'n 
with rriT — m as opposed to px may be a consequence of 
conserving energy instead of momentum in this process, 
though again there is no a priori reason for this to be the 
case. 



fm/c), then the active thermodynamic degrees of free- 
dom during this time period cannot be associated with 
speciflc quasi-particle excitations of the medium. Thus, 
the determination of rj/s of the medium during this early 
time stage is critical to determine whether the medium 
cannot be described as composed of quasi-particles. If 
there are no quasi-particles, then the original concept of 
the quark-gluon plasma as being composed of quarks and 
gluons with well-deflned thermal masses is not realized. 

By comparing the empirical observation of "n^" scal- 
ing of elliptic flow with hydrodynamic fluid calculations 
followed by Cooper-Frye hadronization we come to the 
novel conclusion that at all {rriT — m) ^'riq" scaling gives 
better agreement to the experimental data. It has been 
pointed out that in the region where V2 is linearly related 
to {mT — m), all values of a would preserve v^/a as linear 
with {rriT — rfi)loL. However, as shown in Figure [H typ- 
ical hydrodynamic models with Cooper-Frye freeze-out 
do not evince this linear relationship. 

We also present results for flow scaling for light nu- 
clei and use this comparison to highlight the similarities 
and differences of nucleon coalescence into fragments and 
quark coalescence into hadrons. Just as the light nuclei 
data does not prove that nuclcons were the degrees of 
freedom when the flow developed, the "rig" scaling is not 
sufficient evidence that constituent quarks were the de- 
grees of freedom when the flow developed. 

We find that if one had knowledge of the relevant dy- 
namical masses and quantum numbers of the ensemble of 
quasi-particles it would be possible to formulate a multi- 
stage approach consisting of: hydrodynamic collective 
flow, followed by quasi-particle transport (QPT), coa- 
lescence into hadrons, hadronic re-scattering and flnally 
freeze-out into free streaming hadrons. This proposal of 
a QPT stage between the fluid phase and hadronization 
has yet to be theoretically explored. Most importantly, 
the masses and widths of the quasi-particles will have to 
be related to the creation and total duration of the QPT 
phase. 
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SUMMARY 

In summary, it has been shown that if the matter pro- 
duced in heavy ion collisions is a near-perfect fluid dur- 
ing the stage when the elliptic flow is built up (r < 5 — 7 
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